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i | Abstract. CAT(O) metric spaces constitute a far-reaching common generalization of Euclidean and hyper- 

bolic spaces and simple polygons: any two points x and y of a CAT(O) metric space are connected by a unique 
shortest path ^(x, y). In this paper, we present an efficient algorithm for answering two-point distance queries 
in CAT(O) rectangular complexes and two of theirs subclasses, ramified rectilinear polygons (CAT(O) rectan- 

'— 1 gular complexes in which the links of all vertices are bipartite graphs) and squaregraphs (CAT(O) rectangular 

complexes arising from plane quadrangulations in which all inner vertices have degrees > 4). Namely, we 
show that for a CAT(O) rectangular complex AC with n vertices, one can construct a data structure T> of size 

\{-^ 0(n 2 ) so that, given any two points x, y € AC, the shortest path j(x,y) between x and y can be computed 

in 0(d(p,q)) time, where p and q are vertices of two faces of AC containing the points x and y, respectively, 
such that j(x,y) C K.(I(p,q)) and d(p,q) is the distance between p and q in the underlying graph of AC. If AC 

o 

I is a ramified rectilinear polygon, then one can construct a data structure V of optimal size 0(n) and answer 

two-point shortest path queries in 0(d(p, q) log A) time, where A is the maximal degree of a vertex of G(AC). 
• • Finally, if AC is a squaregraph, then one can construct a data structure V of size 0(n log n) and answer two-point 

• — shortest path queries in 0(d(p,q)) time. 

X 

Keywords. Shortest path problem, rectangular complex, geodesic /2-distance, global nonposi- 
tive curvature. 



1 Introduction 

The shortest path problem is one of the best-known algorithmic problems with many ap- 
plications in routing, robotics, operations research, motion planning, urban transportation, 
and terrain navigation. This fundamental problem was intensively studied both in discrete 
settings like graphs and networks (see, e.g., Ahuja, Magnanti, and Orlin [1]) as well as in 
geometric spaces (simple polygons, polygonal domains with obstacles, polyhedral surfaces, 
terrains; see, e.g., Mitchell [38]). In the case of graphs G = (V,E) in which all edges have 
non- negative lengths, a well-known algorithm of Dijkstra allows us to compute a tree of short- 
est paths from any source vertex to all other vertices of the graph. In simple polygons P 
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endowed with the (intrinsic) geodesic distance, each pair of points p,q G P can be connected 
by a unique shortest path. Several algorithms for computing shortest paths inside a simple 
polygon are known in the literature [29, 30, 33, 36, 43], and all are based on a triangulation 
of P in a preprocessing step (which can be done in linear time due to Chazelle's algorithm 
[12]). The algorithm of Lee and Preparata [36] finds the shortest path between two points of 
a triangulated simple polygon in linear time {two-point shortest path queries). Given a source 
point, the algorithm of Reif and Storer [43] produces in 0(n log n) time a search structure (in 
the form of a shortest path tree) so that the shortest path from any query point to the source 
can be found in time linear in the number of edges of this path (the so-called single-source 
shortest path queries). Guibas et al. [30] return a similar search structure, however their pre- 
processing step takes only linear time once the polygon is triangulated (see Hersberger and 
Snoeyink [32] for a significant simplification of the original algorithm of [30] ) . Finally, Guibas 
and Hersberger [29] showed how to preprocess a triangulated simple polygon P in linear time 
to support shortest-path queries between any two points p, q £ P in time proportional to the 
number of edges of the shortest path between p and q. Note that the last three mentioned 
algorithms also return in O(logn) time the distance between the queried points. In the case 
of shortest path queries in general polygonal domains D with holes, the simplest approach is 
to compute at the preprocessing step the visibility graph of D. Now, given two query points 
p, q, to find a shortest path between p and q in D (this path is no longer unique), it suffices 
to compute this path in the visibility graph of D augmented with two vertices p and q and all 
edges corresponding to vertices of D visible from p or q; for a detailed description of how to 
efficiently construct the visibility graph, see the survey [38] and the book [20] . An alternative 
paradigm is the so-called continuous Dijkstra method, which was first applied to the shortest 
path problem in general polygonal domains by Mitchell [37] and subsequently improved to 
a nearly optimal algorithm by Hershberger and Suri [33]; for an extensive overview of this 
method and related references, see again the survey by Mitchell [38]. 

In this paper, we present an algorithm for efficiently solving two-point shortest path queries 
in CAT(O) rectangular complexes, i.e., rectangular complexes of global non-positive curvature. 
CAT(O) metric spaces have been introduced by M. Gromov in his seminal paper [28] and inves- 
tigated in many recent mathematical papers; in particular, CAT(O) spaces play a vital role in 
geometric group theory [14, 41, 44, 45] (CAT(O) cubical complexes also occur in reconfigurable 
systems [26, 27] and metric graph theory [3]); CAT(O) metric spaces can be characterized as 
the geodesic metric spaces in which any two points can be joined by a unique geodesic short- 
est path, therefore they represent a far-reaching generalization of geodesic metrics in simple 
polygons. Several papers are devoted to algorithmic problems in particular CAT(O) spaces. 
For example, the recent paper by Fletcher et al. [25] investigates algorithmic questions related 
to computing approximate convex hulls and centerpoints of point-sets in the CAT(O) metric 
space P(n) of all positive definite n x n matrices. Billera et al. [8] showed that the space of 
all phylogenetic trees defined on the same set of leaves can be viewed as a CAT(O) cubical 
complex. Subsequently, the question of whether the distance and the shortest path between 
two trees in this CAT(O) space can be computed in polynomial (in the number of leaves) time 
was raised. Recently, Owen and Provan [42] solved this question in the affirmative; the paper 
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[13] reports on the implementation of the algorithm of [42]. 

The remaining part of the paper is organized in the following way. In the next preliminary 
section, we introduce CAT(O) metric spaces, CAT(O) box and rectangular complexes, ramified 
rectilinear polygons, and squaregraphs. We also introduce the two-point shortest path query 
problem. In Section 3, we show that the shortest path 7(3;, y) between two points x, y of a 
CAT(O) box complex JC is contained in the subcomplex induced by the graph interval I(p, q) 
between two vertices p, q belonging to the cells containing x and y, respectively. Moreover, we 
show that this subcomplex )C(I(p, q)) can be unfolded in the /c-dimensional space R k (where k 
is the dimension of a largest cell of JC(I(p, q))) in a such a way that the shortest path between 
any two points is the same in IC(I(p,q)) and in the unfolding of IC(I(p, q)). In Section 4, we 
present the detailed description of the algorithm for answering two-point shortest path queries 
in CAT(O) rectangular complexes and of the data structure V used in this algorithm. First 
we show how to compute the unfolding of )C(I(p,q)) in M 2 efficiently. Then we describe the 
data structure V and show how to use it to compute the boundary paths of fC(I(p, q)). V is 
different for general CAT(O) rectangular complexes, for ramified rectilinear polygons, and for 
squaregraphs. We conclude with a formal description of the algorithm and the analysis of its 
complexity. 

2 Preliminaries 

2.1 CAT(O) metric spaces 

Let (X, d) be a metric space. A geodesic joining two points x and y from X is the image 
of a (continuous) map 7 from a line segment [0, 1] C M to X such that 7(0) = x, j(l) = y 
and d(j(t),j(t')) = \t — t'\ for all t, t' £ [0,/]. The space (X, d) is said to be geodesic if every 
pair of points x,y £ X is joined by a geodesic [11]. A geodesic triangle A(xi,X2,X3) in a 
geodesic metric space (X, d) consists of three distinct points in X (the vertices of A) and a 
geodesic between each pair of vertices (the sides of A). A comparison triangle for A(xi, X2, X3) 
is a triangle A(x / 1 ,x / 2 ,x 3 ) in the Euclidean plane E 2 such that d E 2(x f i ,x'j) = d(xi,Xj) for 
i,j G {1,2,3}. A geodesic metric space (X, d) is defined to be a CAT(0) space [28] if all 
geodesic triangles A(xi,X2,X3) of X satisfy the comparison axiom of Cartan-Alexandrov- 
Toponogov: 

If y is a point on the side of A(xi, x 2 , X3) with vertices x\ and X2 and y' is the unique point on 
the line segment [x^, x' 2 ] of the comparison triangle A(x' l5 x 2 , x 3 ) such that d ¥ 2{x' i , y') = d(xi,y) 
fori = 1,2, then d(x 3 ,y) < d E 2(x' 3 ,y'). 

This simple axiom turns out to be very powerful, because CAT(O) spaces can be characterized 
in several natural ways (for a full account of this theory consult the book [11]). In particular, 
a geodesic metric space (X, d) is CAT(O) if and only if any two points of this space can 
be joined by a unique geodesic. CAT(O) is also equivalent to convexity of the function / : 
[0, 1] — >■ X given by f(t) = d(a(t), f3(t)), for any geodesies a and (3 (which is further equivalent 
to convexity of the neighborhoods of convex sets). This implies that CAT(O) spaces are 
contractible. 
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2.2 CAT(O) rectangular and box complexes 

A rectangular complex JC is a 2-dimensional cell complex JC whose 2-cells are isometric to 
axis-parallel rectangles of the Zi-plane. If all 1-cells of JC have equal length, then we call JC 
a square complex; in this case we may assume without loss of generality that the squares of 
the complex are all unit squares. Square complexes are the 2-dimensional instances of cubical 
complexes, viz. the cell complexes (where cells have finite dimension) in which every cell of 
dimension k is isometric to the unit cube of Analogously, the complexes in which all cells 
are axis-parallel boxes are high-dimensional generalizations of rectangular complexes (we will 
call them box complexes). All complexes occurring in our paper are finite, i.e., they have only 
finitely many cells. 

The O-dimensional faces of a rectangular or box complex JC are called its vertices, forming 
the vertex set V(JC) of JC. The 1-dimensional faces of JC are called the edges of JC, and denoted 
by E(tC). The underlying graph of /C is the graph G(JC) = (V(JC), E{JC)). Conversely, from any 
graph G one can derive a cube (or a box complex) by replacing all subgraphs of G isomorphic 
to cubes of any dimensions by solid cubes (or axis-parallel boxes). We denote any complex 
obtained in this way by ||G|| and call it the geometric realization of G. A cell complex JC is 
called simply connected if it is connected and every continuous mapping of the 1-dimensional 
sphere S 1 into JC can be extended to a continuous mapping of the disk D 2 with boundary 
S 1 into JC. The link of a vertex x in JC is the graph Link(x) whose vertices are the 1-cells 
containing x and where two 1-cells are adjacent if and only if they are contained in a common 
2-cell (see [11] for the notion of link in general polyhedral complexes). Given a subset S of 
vertices of JC, we will denote by JC(S) the subcomplex of JC induced by S and by G(S) (or 
G(JC(S)) the underlying graph of JC(S). 

Computationally, a rectangular complex JC is defined in the following way. Each rectan- 
gular face R of JC is given by the circular list of four vertices and edges incident to R. For 
each vertex v of JC, the neighborhood of v is given as the link graph Link(f); for each edge 
of Link(w) there is a pointer to the unique rectangular face containing v and the edges of 
JC incident to v which define this edge of Link(u). Finally, each point x of JC is given by its 
(local) coordinates in a rectangular cell R(x) of JC containing x (notice that R(x) is unique if 
x belongs to the interior of R(x), otherwise x may belong to several rectangular faces). 

A rectangular or box complex JC can be endowed with several intrinsic metrics [11] trans- 
forming JC into a complete geodesic space. Suppose that inside every cell of JC the distance 
is measured according to an l\- or /2-metric. Then the intrinsic h- or h-metric of JC is de- 
fined by assuming that the distance between two points x, y G JC equals the infimum of the 
lengths of the paths joining them. Here a path in JC from x to y is a sequence P of points 
x = xq, x\. . . x m -i,x m = y such that for each i = 0, . . . , m — 1 there exists a cell Ri of JC 
containing Xi and Xj+i; the length of P is 1{P) = X^S) 1 d(xi, Xj+i), where d{xi,Xi+\) is com- 
puted inside Ri according to the respective metric. We denote the resulting ^i- and ^-metrics 
on JC by d\ and c?2, respectively. 

The interval between two points x, y of a metric space (X, d) is the set I(x, y) = {z £ X : 
d(x, y) = d(x, z) + d(z, y)}; for example, in Euclidean spaces, the interval I(x, y) is the closed 
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line segment having x and y as its endpoints. A subspace Y of a metric space (A, d) is gated 
if for every point x G A there exists a (unique) point x' € Y, the gate of x in 7, such that 
d(x,y) = d(x,x') + d(x',y) for all y £ Y. Gated subspaces are necessarily convex, where a 
subspace Y of (X, d) is called convex if /(x, y) C Y for any x^eF. A half-space H of A is a 
convex subspace with a convex complement. For three points x, y, z of a metric space (A, d), 
let m(x, y, z) = I(x, y) n /(y, 2) D I(z, x). If m(x,y, z) is a singleton for all x, y, z € A, then the 
space A is called median [3, 46] and we usually refer to m(x, y, z) as to the median of x, y, 2 
(here we do not distinguish between the singleton and the corresponding point). A graph G 
is a median graph if (V, dc) is a median space, where da is the standard graph-metric of G. 
Discrete median spaces, in general, can be regarded as median networks: a median network 
is a median graph with weighted edges such that opposite edges in any 4-cycle have the same 
length [2]. Median graphs not containing any induced cube (or cube network, respectively) 
are called cube-free. A median complex is a cube or a box complex of a median graph. We 
will say that a subset S of points of a CAT(O) rectangular complex K, is di-convex for i = 1, 2, 
if S is a convex subset of the metric space (/C, dj) (convex subsets of the underlying graph 
G(JC) will be called graph- convex). 

Now we recall the combinatorial characterization of CAT(O) cubical and box complexes 
given by Gromov. 

Theorem 2.1 [28] A cubical (or a box) polyhedral complex tC with the intrinsic ^-metric di 
is CAT(O) if and only if K, is simply connected and satisfies the following condition: whenever 
three (k + 2) -cubes of K, share a common k-cube and pairwise share common (k + 1) -cubes, 
they are contained in a (k + 3) -cube of fC. 

In some recent papers, CAT(O) cubical polyhedral complexes were called cubings [45]. 
With some abuse of language, we will call cubings all CAT(O) box complexes. The following 
relationship holds between cubings and median polyhedral complexes (this result was used in 
several recent papers in geometric group theory [14, 41]). 

Theorem 2.2 [15, 44] Median complexes and cubings (both equipped with the h-metric) con- 
stitute the same objects. 

In this paper we will mainly investigate the CAT(O) rectangular complexes (i.e., 2- 
dimensional cubings), which can be characterized in the following way: 

Theorem 2.3 [5] For a rectangular complex K, the following conditions are equivalent: 

(i) the underlying graph G(JC) of K, is a cube-free median graph; 

(ii) the metric space ()C,d) is median; 

(Hi) K. equipped with the intrinsic l2-metric c?2 is CAT(O); 

(iv) K. is simply connected and for every vertex x £ V(/C), the graph Link(x) is triangle-free. 
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(c) 

Figure 1: (a) a CAT(O) rectangular complex, (b) a ramified rectilinear polygon, (c) a square- 
graph. 

Typical examples of CAT(O) rectangular complexes (which are illustrated in Fig. 1) are 
the squaregraphs (i.e., the rectangular complexes obtained from the plane graphs in which 
all inner faces are 4-cycles and all inner vertices have degrees > 4 [4, 16]) and the ramified 
rectilinear polygons (i.e., rectangular complexes endowed with the intrinsic /i-metric which 
embeds isometrically into the product of two finite trees [5]). As is established in [5], ramified 
rectilinear polygons are exactly the simply connected rectangular complexes K, in which the 
graph Link(x) is bipartite for each vertex x of /C. 

In median graphs, the halfspaces (the convex sets with convex complements) have a special 
structure and plays an important role. It is well known [7, 39, 40] that median graphs 
isometrically embed into hypercubes. The isometric embedding of a median graph G into a 
(smallest) hypercube coincides with the so-called canonical embedding, which is determined 
by the Djokovic- Winkler relation 6 on the edge set of G : two edges uv and wx are 0- 
related exactly when cIg(u,w) + do{v,x) / dc(u,x) + dc{v,w); see [24, 34]. For a median 
graph this relation is transitive and hence an equivalence relation. It is the transitive closure 
of the "opposite" relation of edges on 4-cycles (i.e., 2-dimensional faces of /C): in fact, any 
two ©-related edges can be connected by a ladder (viz., the Cartesian product of a path 
with K2), and all edges ©-related to some edge uv constitute a cutset @(uv) of the median 
graph, which determines one factor of the canonical hypercube [39]. The cutset Q(uv) defines 
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(a) 



(b) 



Figure 2: Two-point shortest path queries. 

two complementary halfspaces (convex sets with convex complements) W(u, v), W(y, u) of G 
[40, 46], where W(u,v) = {x G X : d(u,x) < d(v,x)} and W(v,u) = V — W(u, v). Conversely, 
for any pair of complementary halfspaces Hi , H2 of a median graph G there exists an edge 
xy such that W(x, y) = H\ and W(y, x) = H2 is the given pair of halfspaces (in fact, all 
edges belonging to the same equivalence 0-class as xy define the same pair of complementary 
halfspaces). 

In this paper, we consider the following shortest path problem in CAT(O) rectangular 
complexes K, endowed with the intrinsic ^-metric d% (we illustrate this formulation in Fig. 2): 

Two-point queries: Given two points x, y of fC, compute the unique shortest path 7(0;, y) 
between x and y in fC. 

3 Geodesies and graph-intervals 

In this section, we show that, given two arbitrary points x, y of a CAT(O) box complex 
JC, the ^-shortest path 7(2:, y) is always contained in the subcomplex induced by the graph 
interval I(p, q) between two vertices p, q belonging to the cells containing x and y, respectively. 
Moreover, we show that this subcomplex IC(I(p, q)) can be unfolded in the ^-dimensional space 
M fc (where k is the dimension of a largest cell of /C(J(p, q))) in a such a way that the ^-shortest 
path between any two points is the same in K,(I(p, q)) and in the unfolding of K,(I(p, q)). 

Proposition 1 If p and q are two vertices of a CAT(O) box complex fC, then fC(I(p,q)) is 
d2-convex and therefore j(x,y) C tC(I(p,q)) for any two points x,y £ K(I(p,q)). 

Proof. According to Theorem 2.2, CAT(O) box complexes are exactly the box complexes 
having median graphs as underlying graphs. Let G = G(fC) be the underlying graph of tC. 
Since G is a median graph, the interval I(p, q) is a convex subset (and therefore a gated subset) 
of G [46]. Additionally, in median graphs each convex set S can be written as an intersection 
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(a) 



(b) 



Figure 3: To the proof of Proposition 1. 

of halfspaces [46] (we will present a simple proof of this fact below - Lemma 3.1). Therefore, it 
suffices to show that for each halfspace H of the graph G, the subcomplex IC(H) is ^-convex. 
Indeed, this will show that fC(S) can be represented as an intersection of ^-convex sets of fC 
and therefore that IC(S) itself is c^-convex. 

Lemma 3.1 A convex set S of a median graph G is the intersection of the half-spaces con- 
taining S. 

Proof. In order to prove the lemma, it suffices to show that for any vertex v of G not 
belonging to S there exists a pair of complementary halfspaces Hi , Hi of G such that v G Hi 
and S C Hi. Since S is convex and G is median, 5 is gated [46]. Let u be the gate of v in 
S (i.e., u G I(v,x) for each vertex x G S). Let v' be a neighbor of v in the interval I(v,u). 
Consider the complementary halfspaces W(v,v') and W(v',v) of G defined by the edge vv' . 
Then obviously v G W(v,v'). On the other hand, since v' G I(v,u) C I(v,x) for any vertex 
a; G 5, by the definition of W(v',v) we conclude that x G W(f/,i>), yielding 5 C W(v',v) and 
concluding the proof. 

Finally, we will prove now that for any pair of complementary halfspaces Hi, H2 of G the 
subcomplexes IC(Hi) and /C(i^2) are ^-convex. Suppose without loss of generality that Hi 
and H2 are defined by the edges of the equivalence class 0j of G. The boundary dHi of Hi 
consists of all ends of edges of 0j belonging to Hi (the boundary BH% of H% is defined in 
a similar way). It was shown by Mulder [39, 40] that in median graphs G the boundaries 
dHi and dHi of complementary halfspaces Hi, Hi induce isomorphic convex, and therefore 
median, subgraphs of G. Hence K,{dHi) and K,(dHi) are isomorphic CAT(O) subcomplexes 
of /C, which we call hyperplanes. Note also that Hi U dHi and Hi U dHi induce convex, and 
therefore median, subgraphs of G. All edges xy of 0j have the same length ij in /C. The CAT(O) 
subcomplex K,{dHi U dHi) of /C is isomorphic to the CAT(O) box complex /Q = (9.ffi x [0, k]. 
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We will show that K{H\) and £(#2) are ^-convex by induction on the number of vertices of 
/C. Suppose by way of contradiction that JC(Hi) is not ofo-convex. Then there exist two points 
x,y G JC(Hi) such that the geodesic 7(2;, y) has a point which does not belong to K,{H{). First 
suppose that 7(2:, y) contains a point z which belongs to KL{H2 \ dHz); see Fig. 3(a) for an 
illustration. Then we can find two points xq G 7(0;, z) D K,{dH.2) and yo £ 7(y 5 -2) l~l K.(dH2). 
Then z G 7(2:0, yo)j showing that K(dH2) is not ^-convex in K, and in the CAT(O) subcomplex 
K.(H2). Since dH2 is a convex subgraph of G and therefore of the median subgraph G(H2) of 
G induced by H2, we conclude that in the CAT(O) complex /C(i?2) not every convex subgraph 
induces a cfe-convex subcomplex, contrary to the induction hypothesis. 

Therefore, we can assume that the geodesic j(x, y) is entirely contained in the subcomplex 
of K. induced by H\\JdH2] this case is illustrated in Fig. 3(b). Since 7(2;, y) is a closed set of /C, 
we can find two (necessarily different) points xq, yo G j(x, y) PI K,(dH\) such that all points of 
the geodesic 7(2:0, yo) except xq and yo all belong to the strip /Q minus the hyperplane K.(dH\). 
Let Cq,C\, . . . ,Ck be the sequence of maximal by inclusion cells of /Q intersected by 7(2:0, 2/0 ) 
labeled in the order in which they intersect 7(2:0, yo) ( so that 2:0 G Co and yo G C&). Let 
Zj G 7(20, yo) H Ci-i n Q, i = 1, . . . , fe, and set zo = a?o, z^+i = yo- Then each of the geodesies 
7(zj_i, Zi) belongs to the cell Q-i, respectively. The intersection Cf of each cell Cj with the 
hyperplane K,(dH\) is a cell of K{dH\) and also a facet of Cj. The orthogonal projection 7Tj 
of each box d (i = 0, . . . , k) on its facet Cf is a non-expansive map (with respect to the &2- 
metric). Notice that vro(2:o) = -20 = yo and 1x^(2,^-1) = Zk+i = yo- On the other hand, since Zi G 
Ci-iDCi and the cells Cj_i and Cj are axis-parallel boxes, we conclude that the two projections 
TVi-i(zi) and 7Tj(zj) are one and the same point z° of the cell Cf_ 1 n C° C Cj_i n Cj. Consider 
the path 7°(xo,yo) between xq and yo in the hyperplane K,(dH\) obtained by concatenating 
the geodesies j(zq = xo, ^i),7(^i, z®), ■ ■ ■ , l{z^ yo = z® +l ). Since for each i = 0, 1, . . . , k, the 
map iTi is non-expansive on Cj, the length of the geodesic j(zf,zf +l ) is less or equal to the 
length of the geodesic j(zi,Zi+±). Therefore 7°(2;o,yo) is a geodesic between xo and yo which 
is completely contained in the hyperplane K,{dH\). Since 7(2:0, yo) H K,{dH\) = {xo,yo}, 
we conclude that xo and yo are connected in /C by two different geodesies 7(2:0, yo) and 
7°(xo,yo), contrary to the assumption that /C is CAT(O). This contradiction establishes that 
the halfspaces of G induce indeed ^-convex subcomplexes of /C, establishing in particular 
that )C(I(p, q)) is c?2-convex for any two vertices p, q of /C. □ 

Proposition 2 If x and y are two arbitrary points of a CAT(O) box complex /C, and R(x) 
and R(y) are two minimal by inclusion cells of K, containing x and y, respectively, then 
7(2:, y) C JC(I(p,q)), where p and q are mutually furthest (in the graph G(K,)) vertices of R(x) 
and R(y). 

Proof. R{x) and R(y) are the unique cells of least dimension such that x belongs to the 
relative interior of R(x) and y belongs to the relative interior of R{y). The sets of vertices of 
R(x) and R(y) are convex, and therefore gated, subsets of G. Let p G R(x) and q G R(y) be two 
mutually furthest vertices of R(x) and R(y), i.e. d(p,q) = max{d(p', q') : p' G V(R(x)),q' G 
V(R(y))}, where all distances d(p', q') are computed according to the graph-distance in G(JC). 
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Since G is bipartite, the choice of the pair p, q implies that all neighbors of p in the graphic 
cube G(R(x)) must be one step closer to q than p, i.e., all these vertices (we denote this set by 
A) belong to the interval I{p,q). Analogously, all neighbors of q in G(R{y)) (we denote this 
set by B) also belong to the interval I(p, q). Since I(p, q) is convex and the convex hull of A in 
G contains the whole graphic cube G(R(x)) while the convex hull of B contains G(R(y)), we 
conclude that both R(x) and R(y) belong to the subcomplex )C(I(p,q)). Since by Proposition 
1, K,{I{p,q)) is a^-convex and x G R(x),y G R(y), we conclude that j(x,y) C )C(I(p,q)). □ 

The next result shows that the intervals I(p, q) in the CAT(O) box complexes can be 
unfolded in Euclidean spaces of dimension equal to the topological dimension (i.e., the least 
dimension of a cell) of JC(I(p,q)). Recall that a function / : X — > X' between two metric 
spaces (X, d) and (X',d') is an isometric embedding of X into X' if d'(f(x),f(y)) = d(x,y) 
for any x, y G X. In this case Y := f{X) is called an (isometric) subspace of X' . If a mapping 
/ : X \— > X' between two geodesic metric spaces (X, d) and (X',d') is such that f{X) is 
geodesic and compact, then f{X) is called an unfolding of X in X' if / is an isometric 
embedding of (X, d) in (f(X),d*), where d* is the intrinsic metric on f{X) induced by d' . 

Proposition 3 If p and q are two vertices of a CAT(O) box complex K, and k is the largest 
dimension of a cell of IC(I(p,q)), then there exists an embedding f* of IC(I(p,q)) in the It- 
dimensional Euclidean space M fc such that f*(JC(I(p,q))) is an unfolding of )C(I(p,q)). 

Proof. First, we show that the subgraph G(I(p,q)) of G = G()C) induced by the interval 
I(p,q) can be isometrically embedded in the /c-dimensional cubical grid Z fc = Il^ =1 Pj, where 
each Pi is the infinite path having 7L as the set of vertices. Indeed, intervals I(p, q) of median 
graphs and median semilattices can be viewed as distributive lattices by setting x A y = 
m(p, x, y) and x V y = m(q, x, y) for any x, y G I(p, q), where m is the median operator of G 
[6, 10]. Using the encoding of distributive lattices via closed subsets of a poset due to Birkhoff 
[9], the famous Dilworth's theorem (the size of a largest antichain of a poset equals to the 
least size of a decomposition of the poset into chains) [21] implies that any distributive lattice 
L of breadth k can be embedded as a sublattice of a product of k chains, see [35] or [19] for 
this interpretation of Dilworth's result (the breadth of a distributive lattice L is equal to the 
largest out- or in-degree of a vertex in the covering graph of L). Larson [35] showed that 
the resulting embedding can be chosen to preserve the covering relation, i.e. to be a graph 
embedding. Recently, using the same tools, Cheng and Suzuki [19] noticed that the embedding 
can be selected to be an isometric embedding of the covering graph of a distributive lattice 
of breadth k in the product of k chains (note that Eppstein [23] showed how to decide in 
polynomial time if a graph G isometrically embeds into the product of k chains). 

Therefore, it remains to show that the largest out-degree or in-degree of G(I(p, q)) equals 
to the dimension of a largest cube of G(I(p, q)). For this, it suffices to show that if v G I(p, q), 
then any m neighbors yi,y2, ■ ■ ■ ,Vm G I(v, q) C I(p : q) of v define an m-cube C rn C I(v, q). 
We proceed by induction on m. Denote by C' the (m — l)-cube induced by the vertices 
yi, . . . , y m -i (which exists because of the induction assumption). Let Zi be the median of the 
triplet yi, y m , q. Then z±, . . . , z m -\ are all adjacent to y m and therefore are pairwise different 
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(because median graphs are .K^-free), and all belong to the interval I(y m ,q). Therefore, 
by induction hypothesis, z±, . . . , z m -i induce an (m — l)-cube C" C I(y m ,q). Then it can 
be easily shown that each vertex of C is adjacent to a unique vertex of C" and that this 
adjacency relation induces an isomorphism between the cubes C and C" . Hence C U C" 
is an m-dimensional cube. Thus indeed G(I(p,q)) can be isometrically embedded in the 
fc-dimensional grid Z fc . 

Denote by / such an embedding. To transform / into an unfolding of JC(I(p, q)) in R k , we 
simply transform the uniform cubical grid Z fc into a non-uniform one: notice that all edges 
of the same equivalence class Gj of G(I(p,q)) are mapped by / to one and the same edge e 
of some path Pj. If the edges of 0j all have length li, then we simply assign length U to the 
edge e of Pj . Denote the resulting paths by Pj* , . . . , P£ and notice that after this scaling the 
previous embedding induce an embedding /* of the graph G(I(p,q)) weighted by the length 
of edges in K, into the grid U^ =1 P*. We can extend in a natural way /* to an embedding of 
K(I(u,v)) into R k = \\U k =1 P*\\ : for a cell R of K(I(u,v)), f*(R) is the cell induced by the 
images under /* of the vertices of R. Let f*(lC(I(p, q))) denote the box complex consisting of 
the images of all cells of IC(I(p, q)). Since each path between two points x,y of IC(I(p,q)) is 
mapped to a path of the same length of f*(IC(I(p,q))) between f*(x) and f*(y), we obtain 
the desired unfolding of K,(I(p,q)) in the /c-dimensional Euclidean space. □ 

For efficient (but nonlinear) algorithms for isometric embeddings of median graphs into 
cubical grids of least dimension, see the recent paper by Cheng [18]. 

4 Two-point shortest path queries 

In this section, we present the detailed description of the algorithm for answering two-point 
shortest path queries in CAT(O) rectangular complexes and of the data structure T> used in 
this algorithm. First we show that IC(I(p, q)) always can be unfolded in the plane as a chain 
of monotone polygons, which we will denote by P(I(p, q)), and we show how to compute this 
unfolding efficiently. Therefore, to compute j(x,y), we triangulate each monotone polygon 
of P(I(p,q)) and compute in linear time the shortest path j*(x,y) in P(I(p,q)) between the 
images of x and y (we denote them also by x, y) using the algorithm of Lee and Preparata 
[36] and return as j(x, y) the preimage of 7*(x, y). As a preprocessing step, we design a data 
structure T> allowing for each query x, y to efficiently retrieve the boundary of an interval 
I(p,q) such that x,y € )C(I(p,q)) and x,p and y, q belongs to common rectangular cells, 
respectively (in time proportional to the distance d(p, q) between p and q in G(JC)). 

4.1 The unfolding of /C(J(p, q)) in R 2 

From Proposition 3, we know that for any two vertices p and q of a CAT(O) rectangular 
complex /C, the graph G(I(p,q)) is isometrically embeddable in Z 2 and consequently the 
subcomplex K.(I(p, q)) can be unfolded in M 2 . We will show how to compute such an unfolding 
efficiently. Denote by P(I(p,q)) the image of IC(I(p,q)) under an unfolding / (which we 
know to exist). Let B\,... ,B m be the 2-connected components (alias blocks) of the graph 
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Figure 4: An unfolding of K(I(p, q)) into blocks. 



G(I(p, q)). From the definition of I(p, q), it follows that B%, . . . , B m define a chain of blocks, 
i.e., if s\ £ B\ n i?2, • • • , %-i S -B m _i fl B m are the articulation vertices of G(I(p, q)) and 
so = p,s m = q, then Bi = J(s$_i, Sj), i = 1, . . . ,m, and all vertices si, . . . , s m _i belong to all 
shortest paths between p and g. In the squaregraph G(I(p,q)) each hyperplane is a convex 
path and in any isometric embedding of G(I(p, q)) in Z 2 , the image of this path is a horizontal 
or a vertical path. Therefore P(I(p,q)) = U^Pj is horizontally and vertically convex (i.e., 
the intersection of P(I(p,q)) with each horizontal or vertical line is convex) and consists 
of a chain of monotone polygons Pi = f(fC(Bi)), . . . ,P m = f(JC(B m )). Up to translations, 
rotations, and symmetries, each block Bi has a unique isometric embedding in Z 2 , whence 
each tC{Bi) has a unique unfolding in the plane. We will define the boundary dG(I(p, q)) 
of G(I(p,q)) (or of K.(I(p,q))) as the subgraph induced by all edges of G(I(p,q)) which are 
mapped to the boundary dP(I(p,q)) of P(I(p,q)) (since G(I(p,q)) is a squaregraph, this 
definition is equivalent to the definition of a boundary of a squaregraph given in [4]). Given 
two arbitrary points x,y £ lC(I(p,q)), it is a well-known property of simple polygons that all 
vertices of the (Euclidean) shortest path in P(I(p,q)) between f(x) and f{y) are vertices of 
the boundary of P( I (p, q)) (except f(x) and f{y)). Therefore, instead of defining the image 
of the whole subcomplex JC(I{p,q)) under an unfolding map /, it suffices to define only the 
image under / of its boundary dG(I(p, q)) (in this case, we will speak about the unfolding of 
dG(I(p,q))). 

Proposition 4 Ifp and q are two vertices of a CAT( 0) rectangular complex JC, then IC(I(p, q)) 
can be unfolded in the Euclidean plane as a chain of monotone polygons and this embedding can 
be constructed in 0(\I(p,q)\) time if I(p,q) is given. Moreover, if the boundary dG(I(p,q)) is 
given together with the degrees dego(z) in G(I(p,q)) of all of its vertices z, then an unfolding 
of dG(I(p,q)) can be constructed in 0(d(p,q)) time, where d(p,q) is the distance between the 
vertices p and q in the graph G(IC). 

Proof. Suppose that either I(p, q) or its boundary dG(I(p, q)) is given and we will show how 
to isometrically embed G(I(p,q)) in Z 2 (the passage from the uniform to the non-uniform 
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grid is the same as in the general case described in the proof of Proposition 3). To embed 
G(I(p,q)), it suffices to detect in linear time the 2-connected components of G(I(p,q)) or 
of its boundary dG(I(p, q)), to embed each block Bi or its boundary dBi, and to compose 
these embeddings into a single one by identifying the images of common articulation vertices 
of blocks (see the discussion preceding this proposition). Therefore we can suppose without 
loss of generality that G(I(p,q)) is 2-connected. Then, up to translations, symmetries, and 
rotations by 90°, G(I(p,q)) has a unique isometric embedding / into the square grid and 
JC(I(p, q)) is embedded in M. 2 as a rectilinear polygon P = P(I(p, q)) which is monotone with 
respect to the coordinate axes. Each vertex of P (including f(p) and /(<?)) is either a convex 
vertex (i.e., the interior angle of P between its two incident edges is equal to 90°) or a reflex 
vertex (i.e., the interior angle between its two incident edges is equal to 270°). Convex vertices 
of P are exactly the images of vertices of degree 2 of G(I(p, q)) (in particular, f(p) and f(q) 
are convex vertices of P), while reflex vertices are images of vertices of degree 4 lying on the 
boundary of G(I(p,q)). Note also that the vertices of degree 3 of dG(I(p,q)) are mapped to 
points lying on sides of P. 

Let 7Ti and 1x2 denote the two disjoint shortest (p, g)-paths constituting the boundary of 
G(I(p, q)). To find the polygon P(I(p, q)) which is the image of fC(I(p, q)) under an isometric 
embedding of G(I(p,q)) in Z 2 , it suffices to find the images of 7Ti and 112 under such an 
embedding. For this, first we scan m and 1T2 in order to detect the convex and the reflex 
vertices of each of them (note that the convex and the reflex vertices on each of these two 
paths alternate). Suppose that convex and reflex vertices subdivide tt\ into the subpaths 
7r{, . . . , -/r™. Then we define the image of p to be the point (0, 0), draw the image of ir\ as a 
vertical path with one end at (0, 0) and having length equal to the length of tt\, then draw the 
image of the path ir 2 as a horizontal path with the beginning at the point where the former 
path ended and of length equal to the length of 7r 2 , and so on, on step i we draw the image of 
the current path ir\ to be orthogonal to the image of the previous path vr^ _1 (the direction on 
ir\ depends on whether vr^ _1 and ir\ share a convex or a reflex vertex). To draw the images 
of the subpaths ir^, ■ ■ ■ , ir™' °f ^ we proceed in the same way but we start by drawing the 
image of horizontal path. This embedding of dG(I(p,q)) extends in a natural way to 

an embedding of G(I(p, q)) : the image of each vertex v lying on the convex path (hyperplane) 
with ends v! € tt\ and it" G -1T2 is a vertex f(u) of the horizontal or vertical path between 
/(it') and f(u") and lying on distance d(v,u') from /(it') and d(v,u") from /(it"). Notice that 
/ is an isometric embedding of G(I(p,q)) into the grid Z 2 because, up to 90°'s rotations, 
there exists a unique isometric embedding in which p is mapped to (0, 0) and this embedding 
necessarily satisfies the properties of / (/ itself is defined in the canonical way). □ 

4.2 The data structure V and the computation of dG{I(p,q)) 

In this subsection, we design the data structure T> for general CAT(0) rectangular complexes, 
for ramified rectilinear polygons, and for squaregraphs, allowing us quickly to compute for two 
arbitrary vertices p, q of /C the boundary paths 7Ti and 1x2 of G(I(p, q)) and the degrees dego(z) 
of all vertices z £ itiU^ = d(G(I(p, q))). The main requirement to V is the trade-off between 
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the space occupied by V and the time for computing m and tt2- Further, we will assume that 
n denotes the number of vertices of /C. Notice that any CAT(O) rectangular complex contains 
0(n) edges and faces. Indeed, \E(K,)\ < 2n and \F()C)\ < n, because any cube-free median 
graph G contains a vertex w of degree at most 2 (if the degree of w is 2, then it belongs 
to a unique rectangular cell of /C, otherwise if the degree of w is 1, then w is incident to a 
unique edge) and, removing w, the resulting graph G' is also cube-free and median. As is 
shown in [5] , as a vertex w of degree at most 2 one can select any furthest vertex from a given 
base-point vertex. Since G' contains n—1 vertices, \E{K,) \ — 1 or \E{K,) \ — 2 edges, and |F(/C)| 
or |.F(/C)| — 1 faces, the required inequalities for G follows by applying induction assumption 
to G' . Finally notice that in all three cases (general CAT(O) rectangular complexes, ramified 
rectilinear polygons, and squaregraphs) at articulation points of G(I(p,q)) the order of paths 
7Ti and 7T2 (which one is left and which one is right looking from p) may change. 

CAT(O) rectangular complexes. In case of general CAT(O) rectangular complexes, using 
Breadth-First-Search (BFS), first we compute the distance matrix D of G(K). Additionally, 
running BFS starting from any vertex u of V(K,) , for each vertex v we compute the list of neigh- 
bors L u (v) of v in the interval I(u,v) (these are exactly the neighbors of v which have been 
labeled by BFS before v). Notice that each list L u (v) contains one or two vertices. Now, T> 
includes the distance matrix D of the underlying graph G(IC) and the lists L u (v), u, v G V{JC). 
T> requires 0(n 2 ) space and can be constructed in 0(|y(/C)| |i?(/C)|) = 0(n 2 ) time. 

Now, we will show how to use V to construct the boundary paths 7Ti and TT2 of G(I(p, q)) in 
0(d(p, q)) time. Before describing the algorithm, first notice that, when requested, the degree 
dego(z) of each vertex z G I(p,q) (and therefore of each z G tt\ U 7^) in the graph G(I(p,q)) 
can be computed in constant time by setting dego(z) := \L p (z)\ + |£q(z)| < 4. To build the 
paths 7Ti and 7T2, we initialize tti := {p} =: iT2- Now suppose that after k > steps, x is the 
last vertex inserted in iri and y is the last vertex inserted in TT2 (if k = 0, then x = p = y). 
Notice that x and y have the same distance k to p and the same distance d(p, q) — k to q. We 
will show how to compute in constant time the neighbor x' of x in tt\ and the neighbor y' of y 
in 7T2 . From the definition of the boundary paths 7Ti and 1T2 follows that each of the edges xx' 
and yy' belongs in JC(I(p,q)) to at most one rectangular cell. To compute x\ we distinguish 
three cases (the vertex y' is computed in a similar way or is defined together with x')\ 

Case 1: x = y. 

If x = y coincides with q, then the algorithm halts and returns the paths n\ and TT2- Otherwise, 
either x = y coincides with p or is an articulation vertex of G(I(p, q)). In both cases, if L q {x) 
consists of a single vertex a, then we set x' := a, y' := a and tt\ := tt\ U {a}, TT2 ■= 1^2 U {o}- 
On the other hand, if L q (x) consists of two different vertices a and b, then set x' := a,y' := b 
and 7ri := tt\ U {a}, tt2 := 1^2 U {b}. 

Now suppose that x / y. 

Case 2: \L q (x)\ = 1, say L q (x) = {a}. 

Then clearly a is the next neighbor of x in tt±, therefore x' := a and tt\ := 7TiU{a}. Additionally, 
if L q (y) = {a}, then set 712 ■= 1^2 U {a} and y' := a (Fig. 5(a)). On the other hand, if 
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(a) deg {b) = 4 (b) deg (b) < 3, b G L q {y) (c) deg (b) < 3, a G L,(y) 

Figure 6: The subcases of Case 3. 



= {a, 6}, then we can set tt2 := it2 U {6} and y' := 6 (Fig. 5(b)). Indeed, since the edge 
ya belongs to two rectangular faces defined by x, a, y, and the median of the triplet p, x, y and 
by y, a, b, and the median of the triplet q, a, b, the vertex a cannot be the neighbor of y in iT2- 
Finally, if L q (y) = {b, c}, then the choice of which of the vertices b and c is y' can be done in 
the same way as the choice of the vertex x' in Case 3 below (Fig. 5(c)). 

Case 3: = 2, say L q (x) = {a, b}. 

First note that since x 1 is one of the vertices a or b, from Lemma 4.1 below follows that the 
case dego(a) = dego(b) = 4 is impossible. Suppose without loss of generality that dego(a) < 3. 
If dego(b) = 4, then Lemma 4.1 implies that x' cannot be the vertex b, thus we can set 
7Ti := 7Ti U {a} and x' := a (Fig. 6(a)) Now suppose that also dego(b) < 3. Then both vertices 
a and b belong to the boundary of K(I(p,q)), one to tti and another one to iT2- Hence, y is 
necessarily adjacent to at least one of the vertices a, b. Notice that y cannot be adjacent to 
both a and b, otherwise the vertices a,b,x,y, and the median of the triplet x,y,p induce a 
forbidden -^2,3 (Fig. 7(a)). Therefore, to decide which of a, b is x' and which is y' it suffices to 
consider to which of the two vertices a, b is adjacent y. If b G L q (y) (and therefore a £ L q (y)), 
then we set n\ := 7riU{a}, x' := a and tt2 := 7T2U{6}, y' := b (Fig. 6(b)). Finally, if a G L q (y) 
(and b £ L q (y)), then we set m := 7Ti U {6}, x' := b and 1^2 '■= ^2 U y' '■= 0, (Fig. 6(c)). 

Lemma 4.1 If \L q (x)\ = 2 and x 1 is the next neighbor of x in tti, then dego(x') < 3. 

Proof. Suppose by way of contradiction that the degree dego(x') of x' in G(I(p,q)) is 4. 
Then necessarily x' is either an articulation vertex or a reflex vertex. We will show that 
in both cases we must have |L g (x)| = 1, contrary to our assumption that |L 9 (x)| = 2. Let 
L q (x) = {a/, x"}. Denote by z the median of the triplet x' , x" , q. Now, if x' is an articulation 
vertex, then necessarily x' must be the closest to q vertex of the 2-connected component 
containing the vertex x. Since the 4-cycle C = (x,x',x", z,x) also belongs to this component 
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(a) L q (x) = {a, b} = L q {y) (b) z = c (c) z = b 



Figure 7: To Case 2 and to the proof of Lemma 4.1. 

and z is closer to q than x' , we obtain a contradiction. Now suppose that x' is a reflex vertex. 
This means that in fC(I(p,q)) x' belongs to three rectangular cells defined by the 4-cycles 
C\ = (x, u, a, x' , x), C2 = (a,v,b,x',a), and C3 = (b,w,c,x' ,6) of G(I(p,q)). The neighbors 
a, x of x' are closer to p than x' and the neighbors b, c of x' are closer to q that x'. Therefore z 
is one of the vertices b or c. In both cases we conclude that the edge xx' belongs in IC(I(p, q)) 
to two rectangular cells C and C\, contrary to the assumption that xx' is an edge of the 
boundary path 7Q (Fig. 7(b) & (c)). □ 

Ramified rectilinear polygons. In case of ramified rectilinear polygons fC, the data structure T> 
consists of an isometric embedding of GQC) into the Cartesian product of two trees Ti, T2 and 
a data structure for nearest common ancestor queries in trees build on each of the tree-factors 
T\ and T^. As was shown in [5], the trees T\ and T2 and the isometric embedding of G := GQC) 
in T\ x T2 can be obtained in the following way. First we construct the equivalence classes 
01, ... , m of the Djokovic- Winkler relation on the edge set of G. As we noticed already, is 
the transitive closure of the "opposite" relation of edges of rectangular cells of /C, and therefore 
the equivalence classes of can be easily constructed in total 0(\F(fC)\ + |-E(/C)| + |V(/C)|) = 
0(n) time. Simultaneously with 0, we can construct the incompatibility graph Inc(G) of G: 
the equivalence classes ©i, . . . , m are the vertices of Inc(G) and two equivalence classes 0j 
and 0j define an edge of Inc(G) if and only if there exists a rectangular face of /C in which 
two opposite edges belong to 0, and two other opposite edges belong to @j (the definition 
of Inc(G) is different but equivalent to that given in [5]). Clearly, Inc(G) can be constructed 
in linear time: it suffices to consider each of the rectangular faces of K, and to define an 
edge between the two equivalence classes sharing edges with this face. The graph Inc(G) is 
bipartite [5], moreover a coloring of vertices of Inc(G) (i.e., of equivalence classes of 0) in 
two colors defines the two tree- factors T\ and Tg. (As is shown in Proposition 1 of [5] this is 
equivalent to a coloring of edges of G in two colors such that opposite edges of each 4-cycle 
C have the same color and incident edges of C have different colors.) 

The trees Tj (i = 1, 2) are obtained from G by collapsing all edges colored with a color 
different from i. Equivalently, to construct Tj, we remove the edges (but keep their ends) 
colored with color i and compute the connected components of the resulting graph G{. Each 
connected component C of Gi defines a vertex of Tj and two connected components C', C" 
of Gi define an edge C'C" of Tj if there exists an edge of G (necessarily colored i) with one 
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Figure 8: Construction of iri and 7T2 in a ramified rectilinear polygon. 



end in C' and another end in C" . Notice that each edge of T% or T2 is labeled by a different 
equivalence class of ©. Now, if u is an arbitrary vertex of G and C\{u) and 6*2(11) are the 
connected components of G\ and G2, respectively, containing the vertex u, then the couple 
f(u) = (Ci(u), 62(11)) are the coordinates of u in the isometric embedding / of G in T\ x T 2 - 
The trees T\ , T2 and the coordinates of the vertices of G in T\ x T2 can be defined in total 
0{n) time. The trees T\ and T2 are preprocessed in linear time to answer in constant time 
lowest common ancestor queries [31]. Additionally, for each vertex v we define the sorted 
list Q(v) of the equivalence classes of to which belong the edges incident to v. These lists 
Q(v),v G V()C), occupy total linear space. These are the three constituents of the data 
structure V. 

Now, we will show how to use T> to construct the boundary paths m and tt2 of G(I(p, q)) 
and the degrees dego(z) in G(I(p, q)) of their vertices z (for an illustration, see Fig. 8). First, 
given the coordinates f(p) = (Ci(p), C2{p)) and f(q) = (Ci(g), Ck(g)) of p and q in T\ x T2, 
using two lowest common ancestor queries, one for C\{jp) and C\{q) in T% and the second one 
for C2(p) and C2(q) in T2, we can compute the path Pi connecting C\{p) and Ci(g) in T\ and 
the path P2 connecting C2{p) and C2{q) in T2 in time proportional to the number of edges on 
these paths. Then G(I(p, q)) can be isometrically embedded in the Cartesian product P\ x P2. 
We will call Pi and P2 vertical and horizontal paths, respectively. We also suppose that f(p) 
and f(q) are respectively the lowest leftmost and the upper rightmost corners of Pi x P2. 

We start by setting 7Ti := {p} =: i\2 and we will construct tti and 1x2 in such a way that 
7Q is the upper path and 1x2 is the lower path of the embedding of G(I(p,q)) in Pi x P 2 . 
Let x be the last vertex of n\ and we will show how to define the next vertex x 1 of Let 
f(x) = (Ci(x),C2(x)). Suppose that the next edge of Pi incident to C\(x) is labeled by 
Qi and the next edge of P2 incident to C2(x) is labeled by @j. Using binary search on the 
sorted list Q{x) we can decide in 0(log(deg(x))) time if x has an incident edge belonging to 
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the equivalence class 0j and/or an incident edge belonging to Qj. If x has an incident edge 
e from 6j, then we set x' to be the end- vertex of e which is different from x (notice that 
C-2,(x') = C-2(x)). In this case, the path m goes vertically. Otherwise, if no edge of Oj is 
incident to x, then necessarily there is an edge e' of Qj incident to x and, in this case, we 
set x' to be the end- vertex of e' which is different from x. Then Ci(x') = C\{x) and the path 
7Ti goes horizontally. Analogously, if y is the last vertex of TT2, in order to define the next 
vertex y' of 1T2, we consider the labels Qi and Qj of the next edges of the paths P\ and P2, 
respectively, and using binary search on Q(y) we decide if y has an incident edge belonging to 
Qj and/or an incident edge belonging to Qi. If y has an incident edge e from Qj, then we set 
y' to be the other end-vertex of e and in this case the path 1x2 goes horizontally. Otherwise, 
if no edge of Qj is incident to y, then there is an edge e' of Qi incident to y and in this case 
we set y' to be the end-vertex of e' different from y and the path P2 goes vertically. 

Notice also that in a similar way we can compute in 0(log(deg(z))) time the degree dego(z) 
in G(I(p,q)) of each vertex z G tt\ U 1T2. As a consequence, the paths 7Ti,7T2 and the degrees 
dego(z), z G 7Ti U 7T2 can be computed in total 0(d(p, q) log A) time, where A is the maximum 
degree of a vertex of G(JC). 

Squaregraphs. Finally, in the case of squaregraphs, as a data structure V we will take an 
encoding of plane graphs defining polygonal complexes of nonpositive curvature presented in 
[17]. Let G = G(JC). The data structure from [17] uses vertex labels of size 0(log 2 n) bits 
(for each vertex of G it uses a label consisting of O(logn) integers of length at most logn) 
and allows for each pair u, v of vertices of G to compute in constant time the distance d(u, v) 
between u and v in G and a neighbor u' of u lying on a shortest path between u and v. 
Additionally, we suppose that the plane graph G is represented as a planar subdivision in the 
form of a doubly-connected edge list [20]. Each edge of G belongs to one or two rectangular 
faces and, using this representation, the vertices and the edges belonging to these faces can 
be listed in constant time. 

Now, given two vertices p and q, in order to construct the boundary paths tt\ and 1x2 of 
G(I(p,q)) and the degrees dego(z) of the vertices z G tt\ U 7T2 in G(I(p,q)) we proceed in the 
following way. The algorithm of [17] returns in 0(1) time a neighbor p' of p in the interval 
I(p,q). Necessarily, p' is a vertex of dG(I(p,q)). Without loss of generality, we assign p' to 
the path m . To find the neighbor p" of p in 112 we consider the edges incident to p in one or 
two faces containing the edge pp' . For each of the end- vertices of these edges which is different 
from p we compute in constant time its distance to q using the algorithm from [17]. If one of 
these vertices is closer to q than p, then we denote it by p" and insert it in 7^. Otherwise, we 
set p" := p' . 

Now suppose that x is the last vertex of ir± and that xq is the vertex preceding x in tt±. 
Suppose also that dego(xo) and the neighbors of xq in G(I(p, q)) have been already computed. 
We will show now how to compute in constant time dego(x), the neighbors of x in G(I(p, q)), 
and the next neighbor x' of x in tt\ (the computation of dego(y), the neighbors of y in 
G(I(p,q)), and the next neighbor y' of y in 712 can be done in a similar way). For this, 
we consider the rectangular faces of K, incident to the edge xqx (there are at most two such 
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faces) and performing in constant time distance queries [17] to q, we compute which vertex 
X\ (if it exists) adjacent to x from these faces belongs to I(p, q) (since the edge xqx belongs 
to 7Ti, at most one such vertex can belong to I(p, q)). Again we distinguish several cases. 

First suppose that the vertex x± exists and x\ 6 I(p,q). Then we consider the face of JC 
(if it exists) incident to xx\ and not containing xq and test if the second neighbor X2 of x 
in this face belongs to I(p,q). If £2 belongs to I(p,q), then again we test if the second face 
incident to the edge XX2 exists and if the neighbor X3 of x in this face belongs to I(p,q). If 
£3 exists and x% G I(p, q), then we set x' := X3 and dego(x) := 4 because in this case we have 
L p (x) = {xo,x\} and L q (x) = {x2,xs}. On the other hand, if X2 exists and X2 £ I(p,q) but 
X3 does not exist or X3 ^ I(p, q), then we set x' := X2 and dego(x) := 3 (in this case, xo,x±,X2 
are the three neighbors of x in I(p,q)). Otherwise, if X2 does not exist or X2 ^ I(p,q), then 
we set dego := 2 and x' := x\ if x\ G L q (x). Finally, if x\ G L p (x), then x is an articulation 
vertex of G(I(p,q)) and we proceed x in the same way as the vertex p. Namely, we set x' 
to be the neighbor of x in I(x,p) returned by the algorithm of [17] (as we noticed already, 
this vertex necessarily belongs to dG(I(p,q))). To compute dego(x) and the neighbors of x 
in I(p,q) we need to find if x has another vertex in L q (x). For this, we consider the edges 
incident to x in the faces incident to the edge xx' . For each of the end- vertices of these edges 
which is different from x' we compute its distance to q using the algorithm from [17]. If one of 
these vertices is closer to q than x, then we include it in L q (x) and set dego(x) := 4, otherwise 
we set dego(x) := 3. 

Now suppose that x\ does not exist or x\ does not belong to I(p,q). Then x is an articu- 
lation vertex of G(I(p, q)). Again, we set x' to be the neighbor of x in I(x, q) returned by the 
algorithm of [17] and check if x has another vertex in L q {x). If such a neighbor exists, then 
we include it in L q (x) and return dego(x) := 3, otherwise we return dego(x) := 2 (in this case, 
L p (x) = {xq} and L q (x) = {x'}). For a given x, each of the operations used to compute x' 
and dego(x) requires constant time, therefore, using the data structure T> of size O(nlogn) we 
can construct the boundary path tt\ and compute the degrees in G(I(p,q)) of its vertices in 
total 0(d(p,q)) time. The path 1x2 can be constructed in a similar way and within the same 
time bounds. 

4.3 The algorithm 

Summarizing the results of the previous subsections, we are ready to present the main steps of 
the algorithm for answering shortest path queries in CAT(O) rectangular complexes, ramified 
rectilinear polygons, and squaregraphs. 
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Algorithm Two-point shortest path queries 

Input: A CAT(O) rectangular complex JC, a data structure V, and two points x, y G JC 

Output: The shortest path 7(2, y) between x and y in JC 

1. Given the rectangular faces containing the points x and y, compute the vertices p, q of JC such that 
x,y e JC{I{p,q)). 

2. Using the data structure V, compute the boundary dG(I(p, q)) and the degrees dego(z) in G(I(p, q)) 
of the vertices z of dG(I(p, q)). 

3. Using the algorithm described in the proof of Proposition 4, compute an unfolding / of dG(I(p, q)). 
Let P(I(p,q)) denote the chain of monotone polygons bounded by f(dG(I(p,q))). 

4. Locate f(x) and f(y) in P(I(p, q)). 

5. Using the algorithm for triangulating monotone polygons (see, for example, [20]), triangulate each 
monotone polygon constituting a block of P(I(p,q)). 

6. In the triangulated polygon P(I(p,q)) run the algorithm of Lee and Preparata [36] and return the 
shortest path J* (f(x), f(y)) = (f(x),z 1 ,...,z m ,f(y)) between f(x) and f(y) in P(I(p,q)), where 
zi, . . . , z m are all vertices of P(I{p, q)). 

7. Return (x, f^ 1 (zi), . . . , f^ 1 (z m ), y) as the shortest path j(x, y) between the points x and y. 



It remains to specify how to implement the steps 1 and 2 of the algorithm. For step 
1, given two points x and y of /C, we are also given two rectangular faces R(x) and R(y) 
containing x and y. Then using the distance matrix D for CAT(O) rectangular complexes, the 
coordinates of the embedding in the case of ramified rectilinear polygons, and the distance 
queries from [17] for squaregraphs, we can compute two furthest vertices p and q, where p 
is a vertex of R(x) and q is a vertex of R(y). This takes constant time because we take the 
maximum of a list of 16 distances between the vertices of R(x) and R(y). In the case of 
ramified rectilinear polygons, to compute the distance d(u, v) in constant time, it suffices in 
each tree Tj (i = 1,2) to keep in V the distance dT t (C,Ri) in Tj from each vertex C to the 
root Ri of Tj. Now, if f(u) = (Ci(u), C^w)) and f(v) = (C\(v), C2(v)), then we compute the 
lowest common ancestor C\ of C\{u) and C\{v) in the tree Ti, the lowest common ancestor C2 
of C2(u) and 62(1?) in T2, and return as d(u,v) the value ((% (Ci(u), R±) + c^i (Ci(v), R±) — 
2d Tl (C 1 ,R 1 )) + (d T2 (C 2 (u), R 2 ) + d T2 (C 2 (v), R 2 ) - 2d T2 (C 2 , R 2 ))- Notice also that the step 4 
of the algorithm also requires constant time: having the coordinates of x in R(x) and of y 
in R(y), since R(x) is the unique rectangular face incident to p in K.(I(p,q)) and i?(y) is the 
unique face incident to q, we can easily locate the images of R(x) and R(y) in the polygon 
P(I(p,q)). Summarizing, here is the main result of this paper: 

Theorem 4.2 Given a CAT(O) rectangular complex K. with n vertices, one can construct a 
data structure V of size 0(n 2 ) so that, given any two points x,y G /C, we can compute the 
shortest l2-path ^(x,y) between x and y in 0(d(p,q)) time, where p and q are vertices of two 
faces ofK. containing the points x and y, respectively, such that j(x, y) C )C(I(p, q)) and d(p, q) 
is the distance between p and q in the graph G(/C). If K, is a ramified rectilinear polygon, then 
one can construct a data structure T> of optimal size 0{n) and answer two-point shortest path 
queries in 0(d(p, q) log A) time, where A is the maximal degree of a vertex of G(K). Finally, 
if K. is a squaregraph, then one can construct a data structure V of size 0(n log n) and answer 
two-point shortest path queries in 0(d{p,q)) time. 
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Open questions: (1) We do not know how to design a subquadratic data structure V allowing 
to perform two-point shortest path queries in CAT(O) rectangular complexes in 0(d(p,q)) 
time or how to use the encoding provided by the isometric embedding of ramified rectilinear 
polygons into products of two trees to remove the logarithmic factor in the query time. 

(2) It will be interesting to generalize our algorithmic results (using Propositions 1-3) to all 
CAT(O) box complexes, in particular to 3-dimensional CAT(O) box complexes. 
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